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1. Introduction 

Lattice QCD is the only ab initio method to compute non-perturbative properties of the bound 
states of quarks and gluons. In this study, we focus on the dynamics of the light degrees of freedom 
inside a heavy-light meson. In the static limit (mg — ► °°), we can properly define and compute 
spatial distributions of matter, of the electric charge and of the axial charge. These can then be 
compared to various quark models which are still the main tool for describing the dynamics of 
excited states. The spatial distributions which we present here constitute a detailed testing ground 
for the quark models, for their choice of parameters, as well as for the basic assumptions, and 
therefore may be a precious means to select among them and/or to improve them. 



2. Static heavy-quark propagator 

The latticized version of the static heavy-quark effective theory Lagrangian (HQET), in the 
frame in which the heavy-quark is at rest, reads [J3j] 

=^HQET = £/i t W \h(x)-\ (x-ofh(x-6)\ + 8m{g 2 )h\x)h{x) 

X 

• The power divergent term 8m(g 2 )h^ (x)h(x) of the discretized HQET is a regularization artifact. 

• Vq(x) is either the lattice gauge link variable in time direction, or a combination of links in 
the neighbourhood of the said link. The former gives the Eichten-Hill action; with the latter, 
depending on the prescription, various actions are obtained, i.e. FAT6, HYP-1, HYP-2, 

The static heavy-quark propagator is given by 

1 + *5(J-v) S (-'" 

n=0 



S h (x,0;y,t) = -^8(x-y)d(t)]JV^n) 



. With the Eichten-Hill action, one immediately realizes that the signal of the two point static- 
light function deteriorates beyond use long before the ground state of the heavy-light meson can be 
reliably extracted. 

FAT6 action was the first proposed to ameriolate the situation. The more recent HYP-1 [|J] 
and HYP-2 [||] actions further reduce the noise to signal ratio (figure |]). In this study, we show that 
HYP- 1/2 can be improved by simply iterating the HYP procedure. 

In short, the HYP procedure consists in decorating a given temporal link by the links situated 
in its neighboring hyper-cube. We refer the reader to the references to HYP- 1/2 above and our 



paper [12] for details. For earlier computations of these distributions see refs. [10]. 



Choice of the static action 



For the spatial distributions we consider, the heavy quark is only a spectator. It stands still 
at a point in space, while we measure the values of our desired local operators at the light-light 
vertex as the light quark is pulled away from the fixed center. This leaves us freedom in the choice 
of static quark action. Judging from figure [j], HYP-2 2 action (HYP-2 iterated twice) is the best. 
Henceforth, results presented are obtained using the that heavy quark action. 
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Figure 1: Plot on the left shows the effective mass obtained using six heavy-quark actions (light quark at K = 0. 1374, 
which is somewhat lighter but close to the physical strange quark mass). The right plot shows the (exponential) growth 
of the statistical error of the two-point function as the time separation between the sources is increased/^/. Notice that 
the error on the correlation function with HYP-2 heavy-quarks grows in the same way as with the HYP-1 2 one. Both 
sink and source are smeared a la Boyle. 

Note also that smearing is highly important for the isolation of ground states in both two/three 
point function calculations. In our study, we adopt the gauge invariant smearing of Boyle For 
the parameters used, we refer the reader to our paper [12]. A more detailed study on the effect of 
smearing will appear in another paper (work in progress). 

3. Distributions 

To extract the desired distributions we computed the following three point correlation functions: 
/ = (1/2) "doublet, i.e. ground states J p = [0~ ,l~],e.g.(B,B*) 



z— fixed 



z.— fixed 



Cy{-t sU t s2 ;r) =^{h(x)ysq{x) q{z + r) T q{z + r) (h{y)y 5 q(y)) 

x,y 

C$ 5 (-t sl ,t s2 ;r) =Y i (h{x)y 5 q{x) q{z, + r) jiYs q(z + ?) (*(y)w(y)) f 

x,y 

j p = (l/2)+doublet, i.e. first excited states J p = [0+, \ + },e.g.(Bl,B{) 

C^\-t sU t s z;r) =£( h(x)lq{x) q{z+r) T q{z + r) (h(y)lq(y)Y 

x,y 

C^l(-t sl ,t s2 ;r) h(x)lq(x) q(z + r) y,-y 5 qiz + r) (ft(y)^75?(y)) t 

The matrix elements are extracted after dividing by the appropriate two point functions 1 

C T {-tsut s2 ;r) 



, r=i,ro 



z.— fixed 



, r=i, 7o 



;— fixed 



Rrr) 



R T (f) 



(3? s ) xexp[-g(t s2 -t si )} 
Cr(-t s i,t s2 ;r) 



3f s ) xexp[-^(4 2 -^i)] 



|f s l|,fs2»0 



(B\qTq(r)\B) = ft^J^r) ,r=l, 7() 

(B\qTq(f)\B*) = f m (r) ,r=yy 5 

(B^qrqir)^) = f % {r) J^r) ,T=l,% 

{Bl\qTq{r)\B { ) = f m (r) S = YiYs 



'We take Zexp(— Et) to be the analytic form of our two point functions. Due to the heavy quark symmetry the 
mesons within each doublet are mass degenerate, and their 2f(£tf)'s are equal too. 
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Figure 2: Radial distribution of the vector, axial and the scalar density within the static-light meson. On the left side 
the distributions are obtained with the external (1/2) ~ mesons, and those on the right with the (l/2) + ones, light-quark 
is K q = 0.1374, which is around the physical strange quark mass. All results are displayed in lattice units. 



ft, fy , and fyys are repectively the matter, electric and axial charge distributions for the (1 /In- 
states. The tilde symbol is added to distinguish the distributions for the (1 /2) + - states (L = 1). The 
results are presented in figure ^ for the case of the light quark being of the mass close to that of the 
physical strange quark. The lattice spacing is about 0. 1 fm. 
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For more details please see ref. [12]. Here we decide to emphasize the following observations: 



• The distributions indeed decrease exponentially when the distance from the static source of color 
grows. 2 By integrating the distributions over spherical volumes, we see that the cooresponding 
charges, as a function of radius, already saturate at around 1 fm. 

poo 

• We checked that the discrete sum Y\ fr(x,y,z) and the integral An \ fl ltted (r) r 2 dr are in 

excellent agreement. This implies that spherical symmetry is well observed. 3 The vanishing of 
the distribution at the boundaries is a necessary condition. 

• Axial and scalar charge distributions of the (l/2) + -doublet change sign between 0.2 and 0.4 fm. 
The other distributions are monotonically decreasing. 

• The vector charge distribution remains of the same form when switching the ground state to the 
excited ones, which is of course consistent with Ward identity. As for the axial and scalar charge 
they are completely different (c.f. fig. ^). 

4. g and g 

g and g are the couplings of the soft pion to the mesons within the same doublet (l/2)~ and 
(l/2) + , respectively. Those couplings are related to the emission/absorption of the P-wave soft 

pion, e.g. B* — > Bn, and B\ — > B°*n. These two axial charges are obtained after integrating the 

r°° , 

distributions presented in the previous section. In particular, g is given by / fyiy$(r) r dr (one 
can use the discrete sum as definition as discussed above.) By using the perturbatively determined 
matching to continuum, our final values are: 

g = 0.526(22) , g = —0.157(48) , linear chiral extrapolation 

gHMChPT = 0.444(18) , gHMChPT = -0.14(6) , chiral perturbation theory §j(4.1) 

where we indicate the two chiral extrapolation procedures. An important remark is that our result 
clearly shows that the axial charges g <C g. This is in contradiction with the claims made on the 
basis of the chiral quark models of refs. ^ that the dynamics within (l/2)~ and (l/2) + doublets 
are very similar This indeed can be inferred from HMChPT too. However, for that to be the 
case the charges g and g should be equal, or at least close to one another, which is clearly not the 
case. 



Notice that g reported here is in excellent agreement with the recent results presented in [ |11[ ] 
(see also references therein). The pionic transition between the heavy-light mesons belonging to 
different doublets is more complicated a problem, which will be addressed elsewhere. 

5. Charge radius of (Pseudo-)scalar mesons 

Electric charge radius is defined by V< r 2 > where, for the ground state meson, 

p oo I poo 

(r 2 } r = J, r 2 Mr)r 2 dr/j o / % (r) r^dr , 



2 Except for /j and fy l y 5 , the logarithm of each distributions can be fitted to a polynomial, /j and fy.y s also die out 
at large r, but to fit them a more elaborate functions are needed, c.f. JI4 I. 
3 Finite lattice spacing only contributes a small correction in this regard. 
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(/ fjo ( r ) r2 dr) 1 is simply the Vector Current renormalization constant Zy obtained non-perturbatively. 
The charge radius for excited state is calculated by substituting fy Q by fy Q . 

In the chiral limit, we obtain (r 2 ) , - = (31.6 ± 0.3) xa 2 ^ (r 2 )""' = 0.334(3) fm 2 . 4 

2 2 

For comparison, the vector meson dominance of the electromagnetic form factor 



(B(p>)\Jtr\B(p))=Ftf)(p + p>) , (/)=6 dFUr) 



dq 2 



(5.1) 



q 2 =0 



with q = p — p' , J s ™ = eQQjnQ+ e q qy^q, and being the electric charge of the heavy(light) quark, 
one obtains that (r 2 ) VMD = 6/m 2 = 0.393 fm 2 . This result is in good agreement with our result. 
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4 We have fixed the lattice spacing to a = 0.0995(4) fm, from ro/a = 4.695(18), as obtained on this same lattice 



and by using rQ hys = 0.467 fm. Charge radius of scalar meson is calculated anogolously. 
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